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Abstract 

We present a simultaneous measurement of the lifetime r^o and mixing parameter Am^. 
We use a sample of about 50000 D*^£^i7£ partially reconstructed decays identified with the 

BABAR detector at the PEP-II storage ring at SLAC, where the flavor of the other B meson is 
determined from the charge of another high momentum lepton in the same event. The preliminary 
results are 

tbo = (1.501 ± 0.008 (stat.) ± 0.030 (syst.)) ps, 
Anid = (0.523 ± 0.004 (stat.) ± 0.007 (syst.)) ps~\ 
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1 INTRODUCTION 



The time evolution of mesons is governed by the overall decay rate T{B^) = and by the 

mass difference Amd of the two mass eigenstates. A precise determination of T{B^) reduces the 
systematic error on the parameters \ Vcb\ and \ Vub\ of the Cabibbo-Kobayashi-Maskawa flavor mixing 
matrix. The parameter ll^dV^^I enters the box diagram that is responsible for B^B^ oscillations and 
can be determined from a measurement of Arrid, with a sizeable systematic error due to theoretical 
uncertainties. 

We describe here a measurement of r^o and Am^ performed using 

^ D*+£~i7i decays^ 

selected from a sample of about 88 million BB events recorded by the BABAR detector at the 
PEP-II asymmetric-energy e~^e~ storage ring, operated at or near the T(4S) resonance. BB pairs 
from the T(4S') decay move along the beam axis with a nominal Lorentz boost {Pj) = 0.55, so that 
the vertices from the two i?-decay points are separated on average by about 260 fim. The B^B^ 
system is produced in a coherent P-wave state, so that flavor oscillation is measurable only relative 
to the decay of the first B meson. Mixed (unmixed) events are selected by the observation of two 
equal (opposite) flavor B meson decays. The probabilities of observing mixed {S~) or unmixed 
(tS"*") events as a function of the proper time difference At are 

-lAtl/r^O 

5^ = (1 ±Pcos(AmdAt)), (1) 

4tbo 

where P is related to the fraction w of events with wrong flavor assignment by the relation D = 

1 — 2w and At is computed from the distance between the two vertices projected along the beam 
direction. 

2 THE BABAR DETECTOR AND DATASET 

We have analyzed a data sample of 81 fb~^ collected by BABAR on the T(4S) resonance, a sample of 
9.6 fb"^ below the resonance, to study the continuum background, and a sample of BB simulated 
events corresponding to about three times the size of the data sample. The simulated events are 
processed through the same analysis chain as the real data. BABAR is a multi-purpose detector, 
described in detail in Ref. pQ. The momentum of charged particles is measured by the tracking 
system, which consists of a silicon vertex tracker (SVT) and a drift chamber (DCH) in a 1.5-T 
magnetic field. The SVT measures the momentum of low transverse momentum charged tracks 
that do not reach the DCH due to bending in the magnetic field. The energy loss in the SVT is used 
in this analysis to discriminate low-momentum pions from electrons. Higher-energy electrons are 
identified from the ratio of the energy of their associated shower in the electromagnetic calorimeter 
(EMC) to their momentum, the transverse profile of the shower, the energy loss in the DCH, and 
the information from the Cherenkov detector (DIRC). The electron identification efficiency within 
the tracking volume is about 90%, and the hadron misidentification probability is less than 1%. 
Muons are identified on the basis of the energy deposit in the EMC and the penetration in the 
instrumented flux return. Muon candidates compatible with the kaon hypothesis in the DIRC are 
rejected. The muon identification efficiency is about 60%, and the hadron misidentification rate is 
about 2%. 

^Charge conjugate states are always implicitly assumed; £ means either electron or muon. 
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3 ANALYSIS METHOD 



3.1 Selection of 5° D*+£-i7i decays 

We select events that have more than four charged tracks. We reduce the contamination from hght- 
quark production in continuum events by requiring the normalized Fox- Wolfram second moment [21 
to be less than 0.5. We select — > D*^£~i7£ events with partial reconstruction of the decay 
D*~^ — > TrfD^, using only the charged lepton and the soft pion (vr^ ) from the D*~^. The is 
not reconstructed, resulting in high reconstruction efficiency. BABAR has already published two 
measurements of r^o |S1 E] and a measurement of sin(2/3 + 7) 5^ based on partial reconstruction. 
This technique was originally applied to B^ — > D*~^£~i7i decays by ARGUS 6 , and then used by 
CLEO 7_, DELPHI 8 , and OPAL [Hj. This is, however, the first simultaneous measurement of 
T^o and Arriii based on partial reconstruction. 

To suppress leptons from several background sources, we use only high momentum leptons, in 
the range 1.3 < p < 2.4 GeV/cJ The vr^ candidates have momenta between 60 and 200 MeV/c. 
Due to the limited phase space available in the D*^ decay, the vr^ is emitted within approximately 
a one-radian full opening angle cone centered about the D*^ flight direction (in the T(45) frame). 
Therefore, we approximate the direction of the D*~^ to be that of the vr^ and estimate the energy 
E£)*+ of the D*~^ as a function of the energy of the vr^ using a third order polynomial, with 
parameters taken from the simulation. We define the square of the missing neutrino mass as 

Ml = i^- Ej,,+ - E,f - (p^.+ + p,)2, 

where we neglect the momentum of the B^ in the T(4S') frame (on average, 0.34 GeV/c), and 
identify its energy with the beam energy ^/s/2 in the e^e~ center-of-mass frame. and are 
the energy and momentum vector of the lepton and y*d*+ is the estimated momentum vector of the 
D*^ . The distribution oi M.^ peaks at zero for signal events, while it is spread over a wide range 
for background events (see Fig. P). 

We determine the B^ decay point from a vertex fit of the £ and tt^ tracks, constrained to the 
beam spot position in the plane perpendicular to the beam axis (the x — y plane). The beam spot 
position and size are determined on a run- by-run basis using two-prong events Its size in the 
horizontal {x) direction is 120 ^m.. Although the beam spot size in the vertical (y) direction is only 
5.6 //m, we use a constraint of 50 //m in the vertex fit to account for the flight distance of the B^ 
in the x — y plane. We reject events for which the probability of the vertex fit, Vy-, does not 
exceed 0.1%. 

We then apply a selection criterion to a combined signal likelihood, calculated from p +, and 
Vvi which results in a signal-to-background ratio of about one-to-one in the signal region defined 
as M-l > —2.5 GeV^/c'*. FigureHshows the distribution of M.1 for the events selected at this stage 
of the analysis. The plot on the top is obtained from the events in which the £ and the vTs have 
opposite charges ( "right-charge" ) , the plot at the bottom from events where £ and the vTs have equal 
charges ("wrong-charge"). We use wrong-charge events as a background control sample, to verify 
that the BB combinatorial background is described well by the simulation. To do this, we add the 
off-peak events, scaled by the ratio of the on-peak to the off-peak luminosities, to the BB Monte 
Carlo scaled to match the data in the region < —4.5 GeV^/c^. We then compare the expected 

^The lepton and pion momenta, the tt^ direction, and E*j^ (see below) are always computed in the T(45') rest 
frame. 
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number of events to the number of data events in the signal region. This ratio is 0.996 ± 0.002, 
consistent with unity. For the rest of the analysis we only consider right-charge events. 



xlO 
1400 



1200 



1000 



800 



600 



BaBar 

preliminary 




M2 (GeV^/c'*) 



Figure 1: Ml distribution for the right-charge (top) and wrong-charge (bottom) events. The 

points correspond to on resonance data. The distributions of continuum events (dark histogram), 
obtained from rescaled off-peak events, and BB combinatorial events (hatched area), obtained from 
the simulation, are overlaid. 



3.2 Tag Vertex and B Flavor Tagging 

We restrict the analysis to events in which another charged lepton ("tagging lepton") is found. To 
reduce contamination from fake leptons and leptons originating from charm decays, we require that 
the momentum of the lepton exceeds 1.0 GeV/c for electrons, 1.1 GeV/c for muons. The tagging 
lepton is used to tag the flavor of the other ("tag") B. The decay point of the tag B is computed 
from the intersection of the tagging lepton track with a beam-spot ellipse centered on the position 
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of the reconstructed B. The beam spot constraint is apphed only in the transverse plane, and 
the width in the vertical direction is inflated to 50/im. We compute the proper time difference 
At between the two vertices from their projected distance along the beam direction (z-axis), Az, 
neglecting the motion in the T(45) rest frame (boost approximation). We obtain At = 
where the boost factor is determined from the continuously measured beam energies. To remove 
badly reconstructed vertices we reject all events with either |Az| > 3 mm or a{Az) > 0.5 mm, 
where a{Az) is the uncertainty on Az, computed for each event. The simulation shows that the 
rms of the difference between the true and measured At is 0.64 ps for 70% of the events, and about 
1.7 ps for the rest. We then select only one right-charge candidate per event according to the 
following procedure: if there is more than one candidate in the event, we choose the one lying in 
the region > —2.5 GeV^/c^. We reject the event if there is more than one candidate (either 
right or wrong charge) in this region. This criterion reduces by about 20% the number of signal 
events. For background studies, we select events in the region < —2.5 GeV^/c^ if there is no 
candidate in the region Ml > -2.5 GeVVc^- We find about 50000 signal events over a background 
of about 27000 events in the data sample. 

3.3 Sample Composition 

Our data sample consists of the following event types, categorized according to their origin and 
to whether or not they exhibit a peak in the A^^ distribution. We consider signal to be any 
combination of a lepton and a charged D* produced in the decay of the same B^ meson. Signal 
consists mainly of B^ D*~^l~V£ decays, with minor contributions from B^ D*^it^£^17£, B^ — > 
D*+T~Vr, B^ D*+D-, B^ D*+DX with r, D^', or D decaying to an and B^ 
D*~^h, with the hadron h misidentified as a muon. A peaking B~ background is due to the 
process B~ D*~^7r~'i~Vi. Non-peaking contributions are due to random combinations of a 
charged lepton candidate and a low-momentum pion candidate, produced either in BB events {BB 
combinatorial) or in e^e~ qq interactions with q = u, d, s or c (continuum). We compute 
the sample composition separately for mixed and unmixed events by fitting the corresponding 
A4l distribution to the sum of four components: continuum, BB combinatorial, B^ — > D*^£~i7i 
decays, and B — > D*TTi~i7£ decays. Due to the production of one or more additional pions, these 
latter events have a different Ai^'^ spectrum from that of the direct process B^ — > D*^i~i7£. We 
measure the continuum contribution from the off-resonance sample, scaled to the luminosity of 
the on-resonance sample. We determine the M-1 distributions for the other event types from the 
simulation, and determine their relative abundance in the selected sample from a fit to the A4l 
distribution for the data. Assuming isospin conservation, we assign two thirds oi B ^ D*itI~i'i 
decays to peaking B~ background and the rest to B^ — > D*TTi~V£ , which we add to the signal. 
We vary this fraction in the study of systematic uncertainties. We assume 50% error on the isospin 
hypothesis that allows us to determine the peaking B~ contribution. 

A possible bias in the Adf, distribution comes from the decay chain B — > £~V£D{X), D Yt:^ , 
where the state Y is so heavy that the charged pion is emitted at low momentum, behaving like a 
TT^ . This possibility has been extensively studied by the CLEO collaboration in Ref. OTJ', where 
the three decay modes most likely to cause this bias have been identified: K* wvr^, K*^ p'^ir^ 
and K*^p^7r^. They quote a systematic error of ±2.3% on the background rate as a result of 
that analysis. Based on their result, we attribute the same systematic error to the number of 
combinatorial events below the signal mass peak. 

Figure 121 shows the fit results for unmixed (left) and mixed (right) events. We use the results 
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of this study to determine the fraction of continuum (/^), BB combinatorial (/^-), and peaking 
B^ background (/^-) as a function of A4f,, separately for mixed (/~) and unmixed (/+) events. 
We parameterize these fractions with polynomial functions, as shown in Fig. |3J 




-10 -8 -6 -4 -2 2 -10 -8 -6 -4 -2 2 

M^2 (GeV^//) M^2 (GeV^//) 



Figure 2: Fit to the Adf, distribution for the unmixed events (left) and mixed events (right). 
includes 5° D*+£~i7i, B° D*+TT°£ue, B° D*+UrT- (r ^ £), B^ D*+D- (D, ^ I), 
and B^ D*^h with the hadron h misidentified as a muon. "i?^"" includes B^ D*^tt^£i'£ and 
B^ D*~7r~^X with the Tr"*" misidentified as a muon. "Other" includes B D*TTVrT (t ^ i), 
B — > D*DX {D £Y), and B — > D*TTh with the hadron h misidentified as a muon. 

3.4 T^o and Arrid Determination 

We fit data and Monte Carlo events with a binned maximum-likelihood method. We divide our 
events into one hundred At bins, spanning the range — 18 ps < At < 18 ps, and twenty a At bins 
between and 3 ps. We assign to all events in each bin the values of At and a/\t corresponding 
to the center of the bin. We fit simultaneously the mixed and unmixed events. We do not use an 
extended likelihood, but instead apply a constraint on the fraction of mixed events, which, for a 
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Figure 3: Fraction of continuum events, peaking B~ , and BB combinatorial background in the 
unmixed (left) and mixed (right) lepton-tagged samples. The fraction of continuum events, 
is parameterized only in the region > —4.5 GeV^/c^. For A^^ < —4.5 GeV^/c^, where just 
continuum and combinatorial backgrounds are present, we assume that /„_ = 0, and we compute 



sample of events with dilution V, reads 



Xd-V+——, (2) 



•'''mix ^•'^ unmix ^ 

where, neglecting the difference AF^ between the two mass eigenstates, the integrated mixing rate 
Xd is related to the product x = Arrid • r^o by the relation 
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We describe the At distribution as the sum of the decay probabihties for signal and background 
events: 



J'^{At,aAt,Ml\TBO,Amd) = f^g{Ml) ■ J'^^iAt, aAt) (4) 

+ f^^{Ml).:F^^{At,aAt) 

where the functions J-f^ represent the probabihty density functions (PDF) for signal (i = B^), 
peaking B~ (i = B~), BB combinatorial (i = BB), and continuum (i = qq) events, and the 
superscript +(— ) applies to unmixed (mixed) events. To account for the ±50% uncertainty on the 
isospin assumption (see Section l3^ . we introduce in the PDF a global scale factor Sb- multiplying 
the fractions /^-, equal for mixed and unmixed events, which we float in the fit but constrained 
to unity with 0.5 error. 

The signal PDF is obtained from Eq. ^ modified to account for the finite resolution of the 
detector. The resolution function is expressed as the sum of three Gaussian functions, defined as 
"narrow", "wide" and "outlier": 

n{5At,aAt) = ^ Vc ' ^^ + ^ e ^^^ + -i^e-^^,(5) 

where 6 At is the difference between the measured and the true value of At, o„ and Ow are offsets, 
the factors Sn and Sw account for possible misestimation of a At- The outlier term, described by 
a Gaussian function of fixed width So and offset Oo, is introduced to describe events with badly 
measured At, and accounts for less than 1% of the events. 

We divide signal events according to the origin of the tag lepton into primary, cascade, and 
decay-side tags. A primary tag is produced in the direct decay B^ I'^VfX. These events are 
described by Eq. ^ with T) close to one (a small deviation from unity is expected due to hadron 
misidentification, leptons from J/V', etc.). We expect small values of o„ and Ow for primary tags, 
because the lepton originates from the B^ decay point. 

Cascade tags, produced in the process B^ —>■ DX, D £Y, are suppressed by the cut on the 
lepton momentum but still exist at a level of 9% as we obtain by floating their relative abundance as 
an additional parameter in the Arud fit on data. The lepton production point is displaced from the 
B^ decay point due to the finite lifetime of charm mesons and the e~^e~ energy asymmetry. This 
results in a significant positive value of the offsets for this category. Compared with the primary 
tag, the cascade lepton usually has the opposite charge correlation with the B^ fiavor. The 
same charge correlation is obtained when the charm meson is produced from the hadronization 
of the virtual W from B^ decay, usually referred to as "Upper Vertex", which results in the 
production of two opposite-flavor charm mesons. We account for these facts by applying Eq. ^ 
to the cascade tag events with negative dilution Va = —(1 — 2fuy) = —0.65 it 0.08, where we 
take from the PDG jllj the ratio of the inclusive semileptonic branching ratios of the b quark 

fuv = ^^fj^ — c^e+)+BR(E c — = 0-1''' 0-04 (the contribution to the dilution from other sources 
associated with the 7rf£ candidate, such as fake hadrons, is negligible). 

Decay-side tags are produced by the semileptonic decay of the unreconstructed D^. Therefore 
they do not carry any information about r^o or Am^. The PDF for both mixed and unmixed 
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contributions is a purely exponential function, with an effective lifetime representing the displace- 
ment of the lepton production point from the decay point due to the finite lifetime of the D^. 
We determine the fraction of these events by fitting the cos 0^+^ distribution (see Fig. ^ , where 
9^+^ is the angle between the soft pion and the tag lepton. Using the results of the cos O^+g fit we 
parameterize the probability for each event to have a decay-side tag as a third order polynomial 
function of cos d +f (see Fig. 
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Figure 4: Distribution of cos 0^+£ for unmixed (left) and mixed (center) events. The points with 
error bars are the data, the shaded histograms show the various sample components after the fit. 
Right plot shows the fraction of tags from decay as a function of cos^ +« for signal events. The 
points with error bars are the data, the solid line is the resulting fitted function and the dashed 
line is the simulation prediction. 

The signal PDF for both mixed and unmixed events consists of the sum of primary, cascade, and 
decay-side tags, each convoluted with its own resolution function. The parameters Sn, S^, So, fw, 
and fo are common to the three terms, while each tag type has different offsets {on,Ow)- To 
reduce the systematic error, some other parameters, in addition to r^o and Am^ are free in the fit, 
namely, all the parameters of the resolution functions, the dilution of the primary tags, the fraction 
of cascade tags, and the effective lifetime of the decay-side tags. We fix the other parameters 
(dilution of cascade tags, fraction of decay-side tags) to their values, obtained as described above, 
and then vary them within their error ranges to assess the corresponding systematic error, as we 
have described above. 

We adopt a similar PDF for peaking B~ background, accounting for primary, cascade, and 
decay terms. Because B~ mesons do not oscillate, we use a pure exponential PDF for the primary 
and cascade tags with lifetime r^- = 1.671ib0.018 ps ^5]. We force the parameters of the resolution 
function to equal those for the corresponding signal term. 
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We describe continuum events with an exponential function convoluted with a three-Gaussian 
resolution function. The mixed and unmixed terms have separate effective lifetime parameters. All 
the parameters of the continuum resolution function are set equal to those of the signal, except for 
the offsets, which arc free in the fit. 

The PDF for combinatorial BB background accounts for oscillating and non-oscillating terms. 
It has the same functional form as the PDF for peaking events, but with independent parameters 
for the oscillation frequency and the lifetimes, which we interpret as effective parameters that do 
not necessarily have a precise physical interpretation. The parameters Sn, S'„,,and /„, are set to 
the same values as those in the signal. 

4 RESULTS 

We first perform the measurement on several Monte Carlo samples. We validate each term of the 
PDF by first fitting signal events, for primary, cascade and decay-side tags separately, and then 
adding them together. We then add peaking B~ background, and finally add the BB combinatorial 
background. We observe the following features: 

• The event selection introduces a bias of (+0.016±0.005) ps on tbo and (-|-0.0027±0.0014) ps"^ 
on Anid- 

• The boost approximation introduces a sizable bias on Am-d (—0.0066 ps~^), as observed by 
fitting the true Az distribution. This bias disappears however when we fit the smeared Az 
and allow for the experimental resolution in the fit function. 

• Introducing B~ peaking background does not introduce any additional bias. 

• Adding combinatorial BB events, we introduce an additional bias of +0.011 ps on t^o and 
of -0.0058 ps on Am^. 

• All the parameters with a clear physical meaning are consistent with their expected value; in 
particular the isospin scale factor Sb- = 1-04 it 0.06 is consistent with unity. 

Based on these observations, we correct the data results by subtracting 0.027 ps from r^o, and 
adding 0.0031 ps~^ to Am^. We assign 100% of the correction as a systematic uncertainty due to 
the observed biases. 

Although we determine the parameters for continuum events directly from the fit to on-peak 
data, we independently fit the off-peak events to verify the consistency with the on-peak continuum 
results. 

We finally perform the fit to the on-peak data. Together with Am^ and r^o , we float most of the 
parameters describing the peaking B~ , BB combinatorial, and continuum background events. The 
values of Am^ and r^o were blinded until completion of the study of the systematic errors. The un- 
blinded fit results are r^o = (1.5280 ±0.0084 (stat.)) ps, and Ama = (0.5200 ±0.0043 (stat.)) ps~\ 
We correct these values for the biases measured in the Monte Carlo simulation, obtaining the 
preliminary BABAR results 

Tbo = (1.501 ±0.008 (stat.) ±0.030 (syst.)) ps, (6) 
Ama = (0.523 ± 0.004 (stat.) ± 0.007 (syst.)) ps~\ 
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The correlation between Am^^ and r^o is —1.2%. Am^^ has sizable correlations with S^- (7.7%) 
and So (7.4%). r^o is correlated with /o (13%) and with a parameter corresponding to the fraction 
of mixed combinatorial B~ events (20%). The complete set of fit parameters is reported in Tabled 

Details on the systematic error are reported in Section |SI Figure El shows the comparison 
between the data and the fit function projected on At, for a sample of events enriched in signal 
by the cut Ml > -2.5 GeVVc"^; Fig. El shows the same comparison for events in the background 
region. Figure [71 shows the plot of the time-dependent asymmetry A{At) — ^unmixedi^t)-NMixed{^t) 
for events in the signal region and events in the background region. 

The agreement between the fitting function and the data distribution is good both in the signal 
and in the background regions. We perform a set of parametrized Monte Carlo experiments to 
assess the quality of the fit. The probability to obtain a lower likelihood is 50%. 
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Figure 5: Distribution of At for unmixed (top) and mixed (bottom) events in the signal 
region with linear (left) and logarithmic (right) scales. The points show the data, the curve is the 
projection of the fit result, and the shaded areas from bottom to top are the contributions from 
continuum, BB combinatorial, peaking B~ with decay-side tag, peaking B~ with primary tag, 
signal with decay-side tag, signal with cascade tag, and signal with primary tag. 



5 SYSTEMATIC UNCERTAINTIES 



The systematic errors are summarized in Table 121 We consider the following sources of systematic 
uncertainty: 

1. Sample composition: the fraction of peaking B~ events in the sample, 5^-, is free in the fit; 
therefore we do not assign any systematic error on this fraction. As mentioned in Section [3.31 
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Table 1: Parameters used in the PDFs from which the hkehhood is calculated. The second column 
shows how they are treated in the fit. The third (fourth) column gives the value employed in the 
data (MC) for the parameters that are fixed or used as a constraint. The last column shows the 
sample in which the parameter is used. The first set of parameters corresponds to peaking events, 
the second to BB combinatorial events, and the third to continuum parameters. The last set refers 
to those parameters of the resolution function that are common to all the data sets. Vl, CI and T>(. 
refer to primary, cascade and decay-side tags, respectively. The symbols o, S, and / correspond 
to offsets, scale factors and fractions in the resolution function; the symbols a and p correspond 
to the fraction of decay-side tags and the fraction of mixed events in the decay-side tag sample, 
respectively. 
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Figure 6: Distribution of At for unmixed (top) and mixed (bottom) events in the background Ad^ 
region with hnear (left) and logarithmic (right) scales. The points show the data, the curve is the 
projection of the fit result, and the shaded areas are the contributions from continuum and BB 
combinatorial background. 

there is a ±2.3% systematic uncertainty on the background rate due to possibly peaking 
combinatorial background. We therefore vary the fraction of BB events in the signal region 
by that amount, repeat the fit, and add the variation in the result to the systematic error 
(entry (a) of Table IS}. We neglect the statistical error on the sample composition because it 
is significantly smaller than this systematic effect. 

2. Analysis bias (entry b): we take 100% of the bias observed in the fit on the Monte Carlo 
sample. 

3. Signal and background PDF description: most of the parameters in the PDF are free in the 
fit and therefore do not contribute to the systematic error. We vary the parameters that are 
fixed in the fit by their uncertainty, repeat the fit, and take the corresponding variation in 
T£o and Anid as systematic errors. We take the uncertainty on r^- (entry (c)), and on Dee 
(entry (d)) from the PDG pT]. 

4. We consider effects due to the detector z scale (entry (e)), the knowledge of the PEP-II boost 
(entry (f)), the actual position of the beam spot (entry (g)), and SVT alignment (entry (h)). 
Detailed studies of these effects have been performed in BABAR in other mixing and lifetime 
analyses (dilepton J2] and fully reconstructed B^ D*~^£~i7£ 14 analyses) and provided 
consistent results. As the methods for vertex reconstruction are very similar we provisionally 
take our systematic error due to these effects from Ref. jl4j . 
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Figure 7: Asymmetry between unmixed and mixed events as a function of At, for events in the 
signal region (left) and in the background region (right). Points with error bars are the data, and 
the curve is a projection of the fit result. 



5. We vary the parameters describing the fraction of decay-side tags by their statistical errors 
(entry (i)). 

6. Binned fitting: we vary the number of bins in At from 100 to 250 and in cxAt from 20 to 50, 
and we repeat the fit. Alternatively, we use the average value of the likelihood in the bin 
instead of the value corresponding to the center of the bin. We take the systematic error to 
be the maximum variation with reference to the default result (entry (j)). 



7. Outlier description: we vary the value of the offset of the outlier Gaussian from 
Alternatively, we use a flat PDF for their description (entry (k)). 



-5 ps to 5 ps. 



8. Fit range: we vary the At fit range from it 18 ps to it 10 ps and the aAt range between 1.8 ps 
and 4.2 ps (entry (1)). 



6 CONSISTENCY CHECKS 

We rely on the assumption that the parameters of the background PDF do not depend on Ai^. We 
verify this assumption for the continuum background with the fit to the off-peak events. To check 
this assumption for the BB combinatorial PDF, we perform several cross checks on the data and 
the Monte Carlo. We compare the simulated BB At distribution in several independent regions of 
Ail with Kolmogorov-Smirnov tests and always obtain a reasonable probability for agreement. We 
fit BB events separately in the signal and background Ml region and compare the parameters of 
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Table 2: Systematic uncertainties. 



Source 


Variation 


6tbo (ps) 


5Amd (ps ^) 


(a) BB fraction 


±2.3% 


±0.001 


±0.001 


(b) Analysis bias 




±0.027 


±0.003 


fc) Tr- 
V / is 


1.671±0.018 ps 


±0.002 


±0.001 


(d) Vrp 


0.65±0.08 


±0.005 


±0.001 


(e) z scale 




±0.006 


±0.002 


(f) PEP-II boost 




±0.002 


±0.001 


(g) Alignment 


- 


±0.006 


±0.003 


(h) Beam spot position 




±0.005 


±0.001 


(i) Decay-side tags 




±0.002 


±0.001 


(j) Binning 




±0.002 


±0.002 


(k) Outlier 




±0.001 


±0.002 


(1) At and (T/\t cut 




±0.008 


±0.003 


Total 




±0.030 


±0.007 



the PDF. We fit the signal plus background Monte Carlo events in the signal region only, fixing all 
the parameters of the BB sample to the values obtained in a fit in the background region, and do 
not see any appreciable deviation from the result of the full fit. Finally, we repeat the fit both on 
the data and the Monte Carlo using different Ai'^ ranges for the background region. Once again, 
we do not observe any appreciable difference in r^o and Am^ relative to the default result. 

7 CONCLUSION 

We have performed a measurement of Am^ and r^o on a sample of 50000 partially reconstructed, 
lepton-tagged B^ D*^i~vi decays. We obtain the following preliminary results: 

tbo = (1.501 ± 0.008 (stat.) ± 0.030 (syst.)) ps, 
Anid = (0.523 ± 0.004 (stat.) ± 0.007 (syst.)) ps"\ 

These preliminary values are consistent with published measurements of r^o and Am^^ performed 
by BABAR with different data sets, and with the world averages computed by the Heavy Flavor 
Averaging Group 12 , as can be seen in table IHl The error we obtain on Am^ is comparable to the 
uncertainty on the present world average. 
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Table 3: Comparison of this result with previous BABAR measurements and with the world average. 
In the case of the world average, we report the total error; in all other cases we show the statistical 
and systematic errors separately. P.R. means partial reconstruction. 



Method 


Arud (ps)~^ 


Tbo (ps) 


Reference 


BABAR Hadronic 


0.516 ±0.016 ±0.010 


1.546 ± 0.032 ± 0.022 


nans] 


BABAR U 


0.493 ±0.012 ±0.009 




[13 


BABAR D*-i+ue 


0.492 ±0.018 ±0.013 


1.523 ±0.024 ±0.022 


m 


BABAR D*'£+ue{P.R.) 




1.529 ±0.012 ±0.029 




BABAR L>*-7r+(P.R.) 




1.533 ± 0.034 ± 0.038 


m 


World Average 


0.502 ± 0.007 


1.536 ±0.014 


pi 


This Measurement 


0.523 ±0.004 ±0.007 


1.501 ±0.008 ±0.030 
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